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Abstract. In 2006 the MAGIC telescope observed 
the well known very high energy (VHE, >80GeV) 
blazars Mrk 421 and Mrk 501 in the course of multi- 
wavelength campaigns, comprising measurements in 
the optical. X-ray and VHE regime. 

MAGIC performed additional snapshot observa- 
tions on Mrk 421 around the MWL campaigns and 
detected the source each night with high significance, 
establishing once more flux variability on nightly 
scales for this object. For certain nights, the integral 
flux exceeded the one of Crab significantly, whereas 
the truly simultaneous observations have been con- 
ducted in a rather low flux state. The MAGIC 
observations contemporaneous to XMM-Newton re- 
vealed clear intra-night variability. No significant 
correlation between the spectral index and the flux 
could be found for the nine days of observations. 

The VHE observations of Mrk 501 have been 
conducted during one of the lowest flux states ever 
measured by MAGIC for this object. The VHE 
and optical light curves do not show significant 
variabiUty, whereas the flux in X-rays increased by 
about 50%. 

In this contribution, the results of the MAGIC 
observations wiU be presented in detail. 
Keywords: VHE MWL blazar 

I. Introduction 

Blazars are radio-loud active galactic nuclei (AGN) 
viewed at small angles between the jet axis and the 
line of sight. Due to relativistic beaming, they are the 
brightest and most variable high energy sources among 
AGNs. BL Lacertae objects, a subclass of blazars |1|, 
are completely dominated by nonthermal emission and 
show only faint or no emission lines at all. Their spectral 
energy distribution (SED) extends over more than 15 
orders of magnitude and exhibits two pronounced peaks 
when plotted as i^F^, versus v. High-peaked BL Lacertae 
objects (HBLs) are characterized by the first peak being 
located in the UV to X-ray range, whereas the second 
peak appears at 7-ray energies (GeV to TeV). 

The non-thermal, beamed emission detected from 
HBLs is commonly ascribed to synchrotron and inverse- 
Compton radiation from ultrarelativistic charged parti- 



cles, accelerated at shocks in jets moving with relativistic 
bulk speed. Progress in understanding the roles of shocks 
and plasma turbulence for particle acceleration and the 
interplay with the dynamics of the jet expansion is 
expected from modelling the SED of blazars. At present, 
it is unclear whether the observed emission is entirely 
due to accelerated electrons, as in Synchrotron-Self- 
Compton (SSC) models (e.g. Q, O, El), or whether 
the high-energy emission is due to pions produced by 
accelerated protons and ions and subsequent pion decay 
(e.g. 0, ISl). 

Owing to their strong spectral and flux variability at 
all wavebands down to minutes (see e.g. 17|), simultane- 
ous multi-wavelength (MWL) campaigns are essential to 
study the SED of HBLs. Of paramount importance are 
hereby the X-ray and 7-ray regimes, where the two SED 
peaks are located for HBLs. In former MWL campaigns, 
the minimum time scale of truly simultaneous observa- 
tions was restricted by the low sensitivity of the par- 
ticipating VHE telescopes, enabling the investigation of 
transient events only in high flux states. With the advent 
of a new generation of imaging atmospheric Cherenkov 
telescopes with improved sensitivity and lowered energy 
threshold such as MAGIC, H.E.S.S. and VERITAS, the 
detection of HBLs in a few hours becomes feasible also 
during low flux states. 

In this proceedings we present MAGIC (ISl, S) 
observations of the well-known TeV blazars Mrk 421 
and MrkSOl in 2006. For individual nights, these were 
conducted (quasi-)simultaneously with H.E.S.S. ifTOl . 
Suzaku ini, XMM-Newton U^, INTEGRAL 113 and 
KVA (14]. 

Mrk 421 (z = 0.030) was the first extragalactic object 
detected at VHE energies ifTsl . Since then, various 
studies from radio wavelength up to VHE 7-rays have 
shown strong variability patterns for this object. Flux 
variations by more than one order of magnitude (e.g. 
1 16 1), and occasional flux doubling times as short as 15 
minutes (|17|, iTsll . IIT9II ) have been reported. Variations 
in the hardness of the TeV 7-ray spectrum during flares 
were reported by several groups (e.g. 1201 . iBTj . IIT6I ). 
Simultaneous observations in the X-ray and VHE bands 
show strong evidence for correlated flux variability ( 1221 . 
II23J, M). 
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TABLE I 

Analysis results for Mrk421 (April and June 2006). 



(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


Observation Night 


S 




T' run 

lime [nj 


-Emin [GeV] 


F{E > Bniin) 


X^rcd, const 


Fq 




X2rod,PL 


2006/04/22 


10.9 


fj 


0.76 


250 


0.917 ±0.112 


1.28/2 


2.37 ±0.33 


2.05 ±0.23 


2.07/2 


2006/04/24 


25.0 


fj 


0.99 


250 


2.32 ±0.13 


2.71/2 


5.14 ±0.40 


2.25 ±0.09 


1.95/3 


2006/04/25 


20.8 


fj 


1.53 


250 


1.34 ±0.09 


1.67/2 


2.99 ±0.31 


2.26 ±0.12 


0.24/3 


2006/04/26 


16.4 


fj 


1.64 


250 


1.08 ±0.09 


1.32/4 


2.37 ±0.26 


2.35 ±0.17 


0.41/2 


2006/04/27 


33.5 


fj 


1.42 


250 


3.21 ±0.15 


1.87/4 


8.05 ±0.42 


2.07 ±0.07 


4.79/4 


2006/04/28 


19.9 


fj 


2.23 


250 


1.14 ±0.08 


4.34/8 


2.38 ±0.22 


2.47 ±0.14 


0.65/2 


2006/04/29 


23.7 


fj 


2.78 


250 


1.04 ±0.06 


41.0/7 


2.35 ±0.17 


2.28 ±0.09 


2.04/4 


2006/04/30 


10.3 


cr 


0.16 


250 


2.39 ±0.33 




6.85 ± 1.06 


1.66 ±0.18 


1.36/1 


2006/06/14 


7.5 


cr 


0.80 


450 


0.341 ± 0.059 


2.36/1 


1.68 ±0.32 


2.38 ±0.44 


1.45/2 



(1) Day at which the observations started, (2) resulting significances S, (3) eltective observation time, (4) minimum spectral energy -fc'min 
used for the calculation of the integral fluxes, (5) integral fluxes F above i?niin (in units of lO^^m^^s"^), (6) fit quality of a constant-flux 
fit to the individual observation nights, (7)-(9) power-law fit results for the differential energy spectra of AN/dE = Fq ■ {E / Eq)~°' with 

Eq = ITeV; Fq in units of IQ-^m-^s-^ 



The second established TeV blazar |24| MrkSOl 
(z = 0.034) has been classified as an 'extreme' HBL, 
having shown a shift of the synchrotron peak of about 
two orders of magnitude to energies higher than lOOkeV 
during observation of a flare by BeppoSAX in 1997 
dHU, |I26]|, ESI). Also at VHE energies surprisingly 
high activity (up to 10 times the Crab Nebula flux) 
has been detected at that time |28|. Reference |^29l 
describes a long term monitoring campaign in 2004, also 
covering the X-ray and TeV energy bands. It confirms 
the presence of a direct correlation between X-ray and 
VHE 7-ray emission, which appears to be stronger when 
the source is brighter. In 2005, when MAGIC observed 
the object between May and July, the source flux varied 
by an order of magnitude. During the two most active 
nights, rapid flux variability with a doubling time of a 
few minutes was observed |7|. Several extensive SED 
studies based on MWL observations of this object were 
reported (e.g. ||3l, ET], lHH, US, ||33l). However, fliese 
observations only covered high flux states, rendering the 
statistics of low state MWL studies sparse. 

II. Observations and Data Analysis 

All MAGIC data were taken during dark nights in 
wobble mode |34|. Observations of Mrk421 have been 
conducted by MAGIC in the nights from MJD 53848 
through 53856 (2006 April 22"^ - 30*'^, accompanied by 
simultaneous KVA observations for most of the nights. 
From MJD 53854.87 - 53855.35 and from 53853.28 - 
53854.27 Suzaku and XMM-Newton, respectively, pro- 
vided additional simultanous X-ray coverage. In the 
Suzaku window, H.E.S.S. conducted observations pre- 
vious to the MAGIC ones, as did Whipple 3.5 h after 
the MAGIC observations in the XMM-Newton time slot. 
Apart from these dedicated campaigns, MAGIC and 
KVA took data on the source at MJD 53901, foflowing 
a targe t-of-opportunity alert from INTEGRAL. Alto- 
gether, the MAGIC pointings resulted in 1 1 .9 h effective 
observation time in April with a zenith distance (zd) 
range of 9° -41° and 0.80h for 43° < zd < 52° in June. 
The detection significances, calculated using formula 17 
of ll35l . were 64.8 cr and 7.5 cr for the April and June 
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Fig. 1. MAGIC LC (E > 250 GeV) for the Mrk421 observations in 
April. The Crab Nebula flux as measured by MAGIC is shown as a 
blue solid line. The dashed red line corresponds to the average flux 
during the low states. Also shown are the observations windows of 
Suzaku (green bar) and XMM-Newton (red bar). 



dataset, respectively (see Table |l]l. For more details see 

ESI- 

The MWL campaign on Mrk501 started on the 
18*'^ of July 2006 (MJD 53935). The nightly MAGIC, 
H.E.S.S. and KVA observations (MJD 53935 - 53937) 
where covered by a continuous Suzaku pointing from 
MJD 53934.77 - 53935.73. The 9.1 h of MAGIC effec- 
tive on-time (11° < zd < 35°) for all three nights resulted 
in a detection of Mrk501 with 13.4 cr significance. 
Observation details can be found in ll37l . 

The data were processed using the analysis and recon- 
struction software package for MAGIC data (EH, ED)- 
Detailed information on the different analysis steps and 
performance are given in |38|. The MAGIC observations 
on Mrk421 presented here are among the first data taken 
after major hardware updates in April 2006, causing a 
higher than usual energy threshold (see ESI)- 

III. Results 

MAGIC detected the sources each night with high 
significance. Independent analyses confirm the results 
presented here. 
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Fig. 2. Differential energy spectra times for the combined low 
(blue) and high (red) state nights of Mrk421. The data were fitted 
using a log-parabolic power law (solid line) and a power law with an 
exponential cut-off (dashed line). 



Fig. 3. Mrk421 intra-night MAGIC LC (E > 250 GeV) for the night 
of the 29"^ of April. The Crab Nebula flux as measured by MAGIC 
is shown as a blue line. 



A. Mrk421 

The results of the daily data analysis are summarized 
in Table I] For Mrk421, not all resulting spectra for the 
individual nights were consistent with a simple power 
law. For MJD 53853 a likelihood ratio test prefered 
a log-parabolic power law 1401 or a power law with 
an exponential cut-off over a simple power law with 
a probability of 94.9% and 95.2%, respectively. Also 
combining all data from April, the likelihood ratio 
probability clearly showed evidence of a parabolic or 
cut-off shape of the spectrum. 

The nocturnal light curve (LC) for the April mea- 
surements is shown in Figure [T] For three nights 
the integral flux of Mrk421 was higher than the 
one of the Crab Nebula, significantly exceeding the 
flux level of the five other nights. The average 
flux in the high state nights (deduced from a sig- 
nal of 46.4 ct) was (2.75 ±0.09) x lO-^m^^g-i, more 
than twice the average flux of the low state nights 
((L09±0.03) X lO-^m-^s-i from 44.6ct). The com- 
bined high and low flux spectra are shown in Figure 
|2] Both could not be described sufficiently by a simple 
power law, but were modelled equally well by a log- 
parabolic power law and a power law with exponen- 
tial cut-off, respectively, yielding frequencies for the 
position of the second SED peak of (297 ±176) GeV 
and (627 ±187) GeV (log-parabolic) and < 250 GeV 
and (603 ±850) GeV (cut-off) for the low and high 
state, respectively. Additionally, the latter model allowed 
us to derive the cut-off energy of the spectrum, lo- 
cated at (3.1 ± 1.5)TeV for the low state nights and at 
(6.3 ±4.2)TeV for the high state nights. 

Thanks to the high significance detections each night 
we were able to derive the VHE LCs for the individ- 
ual nights on sub-hour scale. MAGIC found a clear 
indication of intra-night variability for the observations 
simultaneous to the XMM-Newton pointings, as shown 
in Figure [3] A fit to a constant yielded an unacceptable 
probability of 8.0 x 10"^ %. The LCs of the other nights 




flux at 0.5 TeV [TeV' rrr^ s'] 

Fig. 4. Measured flux at 500 GeV versus spectral index for all 
Mrk421 measurements of this campaign. The blue solid line shows a 
fit of a constant to the data, the dashed red one a linear fit. 

were consistent with a constant flux (see f36l). 

In Figure |4] the flux at 500 GeV is plotted versus the 
spectral index for all nine nights of MAGIC observations 
of Mrk42L A significant correlation of the flux with the 
spectral index is not apparent but cannot be ruled out 
either (x2/ do f const = 17.3/8, x'^/dofunear = 11.5/7). 

B. MrkSOl 

The average spectrum of Mrk501 measured by 
MAGIC in this campaign could be well described by 
a simple power law from 80 GeV to 2 TeV with a 
photon index of 2.79 ± 0.12. The integral flux above 
200GeV was (4.6 ± 0.4) x lO"'^ m^^ ^-i^ correspond- 
ing to 23 % of the Crab Nebula flux as measured by 
MAGIC. 

The light curves obtained from the MWL campaign 
by MAGIC, Suzaku and KVA are shown in Figure |5] 
The measured flux in the VHE regime did not vary 
significantly and was comparable with the lowest value 
of MAGIC measurements in 2005 and 2006. The Suzaku 
XIS data showed a flux increase of about 50% during 
the observations. 
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Fig. 5. MAGIC (top), Suzaku (middle) and KVA (bottom) LC 
for Mrk501 derived in thi.s campaign. The dashed horizontal lines 
represent the average fluxes of the measurements. 



IV. Discussion 

The well-know TeV blazai's Mrk421 and MrkSOl 
have been observed by MAGIC in the course of several 
MWL campaigns. Most of the simultaneous observations 
were prearranged and catched the sources in rather low 
flux states, adding statistics to the small sample of truly 
contemporaneous low state data. 

Whereas the monitoring of Mrk421 througout the 
campaigns exhibited strong variability between the in- 
dividual observation nights as well as intra-night, the 
observations on MrkSOl revealed no significant varabil- 
ity. The intra-night variability pattern found for Mrk421 
did not show up in the contemporaneous X-ray data by 
XMM-Newton (see Ell). On the other hand, the flux 
increase of Mrk 50 1 in the soft X-rays could not be found 
in the VHE data. However, in this second case, due to 
the low source flux level, MAGIC could only have seen 
variability if the flux were to increase by a factor of 2-3, 
a level of change which was also not apparent in X-rays. 

Using the results of the nine days of MAGIC obser- 
vations on Mrk 421 we searched for correlation between 
the fluxes and spectral indices. No significant correlation 
was apparent, which may be due to the small dynamical 
range of the parameters found here. 

Dividing the observations in a low and high state data 
set, we found that both spectra were not compatible 
with a simple power law but were well fit by a log- 
parabolic power law and also a power law with expo- 
nential cut-off. Our data did not allow to prefer one 
model over the other If source-intrinsic, a log-parabolic 
spectral shape in the VHE regime can be explained by 
stochastic, energy-dependent acceleration of electrons 
within the SSC framework ll40l . Electron scattering in 
the Klein-Nishina regime or absorption of the 7-rays by 
extragalactic background photon fields (e.g. [42]) can 
be responsible for an intrinsic or cosmological spectral 
cut-off feature, respectively. 

The peak frequencies dervied from both fits were 
compatible with each other and showed indications for 
an increase of the peak energy with rising flux level. 



as predicted if the VHE radiation were due to SSC 
mechanisms. The derived cut-off energies were well 
below the limit of ~ 13 TeV deduced from the Kneiske 
lower limit model |42|, indicating a source-intrinsic 
rather than a cosmological reason for the cut-off feature. 

A detailed discussion of the data in a MWL context 
will be presented in subsequent publications, but some 
aspects can already be found in 1411 and i37i . 
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